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Abstract: In this article we discuss a general effective-theory description of a multi- 

component dark sector with an unspecified non-trivial symmetry and its interactions with 
the Standard Model generated by the exchange of heavy mediators. We then categorize the 
relevant effective operators given the current experimental sensistivity where the underly¬ 
ing theory is weakly coupled and renormalizable, with neutral mediators: either scalars or 
fermions. An interesting scenario resulting from this categorization is the neutrino portal, 
where only fermion mediators are present, and where the dark sector consists of fermions 
and scalars such that the lightest dark particle is a fermion, this scenario is characterized by 
having naturally suppressed couplings of the DM to all SM particles except the neutrinos and 
has received little attention in the literature. We find that there is a wide region in parameter 
space allowed by the current experimental constraints (relic abundance, direct and indirect 
detection limits); the cleanest signature of this paradigm is the presence of monochromatic 
neutrino lines with energy equal to that of the DM mass, but experimental sensitivity would 
have to be improved significantly before this can be probed. 


Keywords: Effective Theory, Dark Matter, Neutrino physics 



Contents 


1 Introduction 1 

2 DM-SM interactions 2 

2.1 Tree-level and loop-level generated operators - neutral mediators 4 

2.2 The effective Lagrangian 5 

2.2.1 Some simplifications 6 

3 Scalar mediators 7 

4 Fermion mediators 7 

5 Fermion DM with fermion mediators: neutrino portal DM 8 

5.1 Relic abundance 8 

5.2 Direct detection 11 

5.3 Indirect detection 13 

5.4 Collider onstraints 14 

6 Conclusions 16 

A The operators 17 

B Cross sections 17 


1 Introduction 

Understanding the fundamental nature of Dark Matter (DM) is one of the most compelling 
problems in particle physics and cosmology [1], yet despite significant and continuing exper¬ 
imental efforts, no information about the nature of DM has been obtained (excepting its 
gravitational effects [2-5]). In the current favored paradigm DM is composed of one or more 
particles whose interactions with the Standard Model (SM) are weak enough to meet the 
constraints of direct and indirect detection experiments [6-9], but strong enough to generate 
the relic abundance inferred from measurements of the cosmic background radiation [10] . 

The absence of information on the interactions of the dark sector with the SM, indicates 
that a model-independent study of these interactions will be useful in understanding the 
effects of the various possible couplings, and of the processes that generate them; (see, for 
example, [11-13]). We will follow such an approach by using an effective Lagrangian to 
parametrize the interactions of the dark sector with the SM, and determine the restrictions 
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imposed by the above constraints. We will consider a general dark sector that can include 
vector, Dirac fermions and scalar particles, and these interact with the SM particles through 
the exchange of heavy mediators that we denote generically by 3 ^ . We will assume that the 
mediators are weakly coupled to both the standard and dark sectors and that they satisfy 
the requirements of the decoupling theorem [14], in particular the interactions they generate 
between the dark and standard sectors vanish as the mediator mass —?■ 00 . 

We will assume that all dark fields transform non-trivially under a symmetry group ^dm 
( whose nature we will not need to specify), while all SM particles are assumed to be ^dm 
singlets; these characteristics provide a simple way to ensure the dark sector contains a stable 
particle that will play the role of DM. Finally, we assume that all dark fields are singlets 
under the SM gauge group ^SM = SU{3) x SU{2) x U{1). 

This paper is organized as follows. In the following section we construct the leading 
terms in the effective interactions between the standard and dark sectors, we discuss briefly 
the hierarchy of these couplings generated by the canonical dimension of the corresponding 
effective operator and by whether it is generated at tree-level or not; we use this to identify 
the most phenomenologically interesting interactions. We then specialize to the case where 
the mediators are scalars or fermions and are singlets under ^dm x ^sMj which we consider 
separately in sections 3 and 4. We then concentrate on a specific scenario, which has received 
limited attention in the literature: that of a multicomponent dark sector, with fermionic DM 
and fermionic mediators (section 5). Parting comments are provided in section 6 while details 
of the calculations are relegated to the appendices. 

2 DM-SM interactions 

An immediate consequence of our assuming mediator-generated interactions between the 
standard and dark sector is that these interactinos take the form 

O = OsMOdark (2.1) 

where OsM,dark denote local operators composed of standard and dark fields respectively, 
and which are invariant under both ^dm and ^smj but need not be Lorentz singlets. The 
operators O will apepar in the effective Lagrangian multiplied by coefficients proportional 
to l/Mg with n = dim(O) — 4; in particular, the larger the dimension of O, the smaller its 
effect, so that within this paradigm, dominant effects will be generated by lower dimensional 
operators. Given the detailed knowledge of the SM constructing the operators Osm invariant 
under ^sm x ^dm is a straightforward exercise 

^In this context, ‘heavy’ indicates that the mediator masses Mg are assumed much larger than the typical 
energies of processes involving interactions between the standard and dark sectors. The standard and dark 
sectors presumably also interact via the exchange of gravitons; we assume the effects of such interactions are 
small. 

^Invariance under Qdm is automatic since all SM fields are assumed to be Qdm singlets. 
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VII 

TsM Vdark 


VIII 

'°/^i"'^dark 


Table 1. List of operators of dimension < 6 involving dark and SM fields; (p denotes the SM scalar 
isodoublet, B the hypercharge gauge field, and £ a left-handed lepton isodoublet; also, e = ftT 2 , where 
tT 2 denotes the usual Pauli matrix. Dark scalars, Dirac dark fermions and vectors are denoted by $, 
and X respectively. The currents in category VII operators are defined in (2.2), and the operators 

0 ( 4 ) 

in categories V and VIII are listed in appendix A. See the text for details. 


When constructing Odark we will assume that the dark sector is composed of scalars <k, 
Dirac fermions 'k and vectors JV, with the understanding that the dark sector present in 
Nature may only contain a subset of these particles. Again by assumption all dark fields 
are invariant under ^sMj so that Odark will also be a singlet under this group. Whether a 
given combination of dark fields is invariant under ^dm or not depends on the details of this 
symmetry and the representations carried by these fields, unfortunately existing data does 
not provide any information on this point, so we opt for the most general case and study the 
effects of all interactions of the form ( 2 . 1 ) when Odark is invariant under at least one choice of 
^dm; we can only say that since all dark fields are assumed to transform non-trivially under 
^DM, Odark will Contain at least two fields. The list of the operators O of dimension < 6 that 
satisfy the above conditions is given in table 1 . 

Referring to this list we make the following observations: 

• Some of the operators might be disallowed by the spectrum in the dark sector (e.g. the 
operator in category I would be absent if there are no dark scalars). Other operators 
are present only when the dark sector has several components (e.g. dark fermions and 
scalars are required for the operator in category III to be present). 

• Though it is customary to assume ^dm is a discrete symmetry (e.g. Z^,) this need not be 
the case. For example, ^dm could be a gauge group, and X the associated non-Abelian 
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(dark) gauge bosons. In this case operators such as would be invariant if the 

corresponding belong to the adjoint representation of ^dm- 

• Some operators might be disallowed by the choice of ^dm; for example, if ^dm = ^ 2 , 

our requirement that all fields transform non-trivially under this symmetry requires 
them all to be odd under Z 2 , hence is forbidden. 

• The in categories V and VIII represent dimension-4 field combinations of the 
corresponding sector, invariant under ^dm and ^sMi they are listed in appendix A. 

• In category VII, JsM.dark represent SM and dark currents of dimension 3, invariant 
under ^dm and ^sm: 


AT" = 4*1" = A-I-* A $ 

where '0 denotes any SM fermion, D the covariant derivative in the standard sector, 
and T) the covariant derivative in the dark sector (replaced by an ordinary derivative if 
this sector is not gauged). 

• In general there can be more than one dark field of each type, but we did not include such 
‘dark-flavor’ indices for notational simplicity; similarly we did not include a generation 
index in the left-handed lepton isodoublet i (categories III and VI). 

• There are no operators of the form because they would not invariant under 

^DM) given our assumption that all dark fields transform non trivially under this group. 

2.1 Tree-level and loop-level generated operators — neutral mediators 

All the operators in table 1 will, in general, be generated in any model within the class of 
theories being considered here. Depending on the details of the model, any given operator will 
be generated at tree level, or at one (or higher loops) involving a combination of mediator, dark 
and standard internal lines Operators that are generated at L-loops will have coefficients 
suppressed ^ by a factor ~ l/(47r)^^; tree-generated operators will have no such suppression. 
Because of this, evaluating the effects of an operator and eliciting constraints from data on 
the model parameters will depend strongly on whether the operator is tree or loop generated, 
and this requires that some properties of the mediator sector be specified. 

We will now restrict ourselves to models where the dark and a standard sector interact 
through the exchange of neutral mediators, that is, which are singlets under ^dm x ^sm- 

^Our effective theory will be valid at scales below the mediator mass Mg, and is generated by integrating 
out the mediators as well as all standard and dark field modes with momenta above Mg. 

"^The presence of loop-snppression factors does not necessarily mean that the effects of the corresponding 
operators are phenomenologically irrelevant since the experimental sensitivity that may be sufficiently high to 
warrant retaining them; ignoring such factors, however, can lead to significant over-estimation of these effects. 


( 2 . 2 ) 
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Table 2. Dark-SM operators that are generated at tree level by neutral (under ^dm and ^sm) scalar 
(S) and fermion (J^) mediators. Some of the operators in category V (see table 1) are not listed 
because they are generated at tree-level only when heavy charged (under the dark or SM symmetries) 
particles are present. The currents JsM.dark are assumed to be conserved (see text). 

We will also assume that the full theory composed of mediators, dark and standard sector 
is renormalizable and that the mediators are weakly coupled. Within this neutral-mediator 
paradigm one can determine by inspection that the operators in table 1 are generated at tree 
level by scalar mediators S (categories II an V), fermion mediators J- (categories III and VI), 
vector mediators (category VII), or antisymmetric tensor mediators (categories IV and VIII). 
In this paper we will consider only the case of scalar and fermion mediators; tensor and vector 
mediators require a discussion of the corresponding local symmetries and will be presented 
elsewhere. The operators that can be generated at tree-level by scalar and fermion neutral 
mediators are listed in table 2 

2.2 The effective Lagrangian 

For processes with typical energies below the mediator mass the relevant physics is described 
by an effective theory resulting from integrating out all modes with energies above this scale; 
aside from the mediators themselves, these modes include those of the standard model and 
dark fields with momenta above M^. In the following we refer to these as ‘high-momentum 
modes’ (HMM). Integrating out the 3 and HMM generates an effective Lagrangian of the 
form 

= ^SM + W + + £(3-tree) ^ ^(3-loop) (2.3) 

where the first and second terms correspond to the Lagrangians for the standard and dark 
sectors (including the loop-generated interactions resulting from integrating out the HMM 
within the corresponding sector). The third term contains the renormalizable Higgs-portal 
interaction(s) in category I, where Ci ~ 1 is the natural value of this coefficient (this term is 
of course absent if the dark sector does not contain scalars). 

The fourth term contains the tree-level generated effective dark-SM interactions resulting 
form integrating out the mediators only; these are listed in table 1 . The last term contains 
all loop-generated terms obtained by integrating out mediators and HMM, they generate 
the remaining terms in table I but the corresponding coefficients suppressed by loop factors 
~ l/( 47 r)^'^, where L is the number of loops. The interactions contained in these last two 
terms will be suppressed by inverse powers of the mediator mass ® scale M^. 

is the same scale used in integrating-out the HMM. 
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VII 

|,/.| 2 Tc 1 >T', 4>2(IT^)2 
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Table 3. Leading dark-SM interactions of dimensions 5 and 6 (see section 2.2.1); the operators in 
categories III and VII should be removed when only scalar mediators are present. 

2.2.1 Some simplifications 

For the case where the SM and dark sectors are connected only by exchange of neutral media¬ 
tors the effects of several operators in table 1 are either very small effects or subdominant. In 
the following we will ignore such contributions to the effective Lagrangian (recognizing that 
in near-future experiments they might generate small but potentially observable deviations 
from the leading effects). The operators that we will drop are 

• Operators of dimension > 5 with only scalar fields. Given the current experimental 
sensitivity, the effects of operators in categories II and V that involve only scalars (</> 
and <h; see also appendix A) will be subdominant compared to those generated by the 
Higgs portal coupling oc Ci in (5.6), provided all coefficients are within their natural 
ranges. 

• Operators in category VI. The observable effects generated by these operators are very 
similar but subdominant to those generated by the operator in category III. 

• Operators generated at >2 loops. Such operators appear multiplied by a small coeffi¬ 
cient ~ l/( 47 r)‘^ ~ 4 X 10“^. Specifically, these are: 

— Category III operators when only S mediators are present. 

— All category IV and VIII operators. 

— The operators |‘hp(G)^^,)^ in category V. 

— Category VII operators when only S mediators are present, or when IF mediators 
occur and the operators involve Jgj(f for i ^ i, cp. 

The remaining operators are listed in table 3. 

®These considerations depend strongly on our assumptions of weak coupling and renormalizability. 
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3 Scalar mediators 


When scalar mediators are present, the list of tree-generated operators is on the first line of 
table 2; of the remaining operators in table 1 we retain only those generated at one loop. The 
resulting effective Lagrangian is 


^(5_treee) ^ + 

£(5-loop) _ g(6) . . . 

^ (47rA)2 ^ 


(3.1) 


where the prime indicates that the sum is over the category V operators listed in table 3. 

In the case where both dark scalars and fermions are present the S will also generate the 
dimension 5 interaction |4>p(TT) so that, only the lightest of these particles will contribute 
significantly to the relic density (unless the mass splitting is small or some coefficients are 
significantly suppressed [15-17]). It then follows that, as far as direct and indirect detection 
are concerned one can consider a single-component model where only the lightest of the 
dark particles is included. Such single-component DM models are obtained by retaining 
the appropriate subset of the above interactions. In particular, the Higgs-portal coupling 
(proportional to Ci) is the most significant interaction in models where only dark scalars are 
present; this type of models have been extensively studied in the literature [18-23] . 


4 Fermion mediators 

When only fermion mediators are present ^ the effective Lagrangian resting form integrating 
the J- takes the form 

tree) 

£(J'-loop) 


where the prime indicates that the sum is over the category V operators listed in table 3, 
and the ellipses denote subdominant operators (see section 2.2.1); the currents are defined 
in (2.2). It is also worth noting that the category II operator has a loop-suppressed 

coefficient for these models, which was not the case when scalar mediators are present. 

When expanded in terms of component helds in the unitary gauge (for the SM) the 
dimension 5 term in (4.1) is seen to contain a vertex of the form Cni7'L4>^'I', where ul denotes 
a SM left-handed neutrino. The presence of this coupling implies that the heavier of the 

^When both scalar and fermion neutral mediators are present the effective theory is obtained by adding 
(3.1) and (4.1). The resulting effective theory has an 0(1) coupling for the operators in categories II and III, 
and 0(l/167r^) for categories V and VII. 




Cii 


lOvr^A 




E 

a=i (j)] i=L^R^^ 
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(4.1) 
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dark particles (‘k, dt) will decay promptly into the lighter one and a neutrino. So this type of 
models, while having a multi-component dark sector, have a single-component DM relic 

If the relic DM is composed solely of scalars ‘k and the model largely reduces 

to the one discussed in the simplest Higgs-portal models (when some of the coefficients are 
suppressed by (dvr)^). 

If, however, < m#, the situation is different: DM will be composed of the dark 
fermions dt, and the leading 'k interactions are with neutrinos or neutrinos and Higgs particles. 
This is a case that we discuss in more detail below, and which has not received much attention 
in the literature (right-handed neutrino portals with lepton number violation were examined, 
for example in [24, 25]; and models including heavy Majorana right-handed neutrino mediators 
and fermion DM relics were proposed to explain dark matter-antimatter asymmetry through 
leptogenesis [24, 26, 27]). 

5 Fermion DM with fermion mediators: neutrino portal DM 

We now consider in some detail the case where the dark sector contains Dirac fermions and 
scalars, with the latter being heavier than the former, and when dark particles interact with 
the SM via the exchange of neutral Dirac fermions conserving lepton number ® . As indicated 
above, the main DM-SM interactions involve neutrinos or Higgs particles, while all other 
interactions have loop-suppressed couplings. 

5.1 Relic abundance 

The leading DM-SM interaction is generated by the (tree-level) exchange of the dark scalars 
<k and represent the most important reaction responsible for the equilibration between the 
dark and standard sectors in the early universe; this process is produced by the interaction 
oc Cm in (4.1). It is worth noting that even if the dark scalars do not contribute to the relic 
abundance, their presence is essential for the viability of the model: in the absence of 4> all 
terms in (4.1) would be absent and the SM and dark sectors would decouple. 

The remaining interactions in £(-^-1°°?) generate small corrections but for two exceptions: 
the terms oc Cn and oc contain the vertices TTR (where H denotes the Higgs) and 

^, and will produce important resonant effects when ~ fni{/2 and ~ Tnzj^ 

that are generated, respectively, by the s-channel exchange of the H and Z. Other interactions 

® Since the neutrino masses are so small there is no need to assume a large splitting between m<j and 
to differentiate between these two scenarios. The more complicated case where the dark fermions and scalars 
are degenerate is unnatural as there is no symmetry that can ensure this constraint. For example, if the dark 
sector is assumed to be supersymmetric with >1> and 'k members of a multiplet, then the category III operator 
explicitly breaks this dark supersymmetry. 

®When model-building, small violations of lepton number can be introduced in the mediator sector as 
perturbations. This would not affect the discussion below in any significant (qualitative and quantitative) 
way. 

^°For the loop-suppressed terms this follows from a straightforward examination of the possible diagrams. 



of the form Jsm • <^dark or ill (4-1) generate small corrections that we will ignore in the 
following. 

In the unitary gauge 

, ^^\cj,\‘^DvH{m)+, 4M-4t;S^ (5-1) 

with V ~ 246GeV. The cross section for —)• i/v (relevant for the relic abundance calculation 
below) is generated by the diagrams in figure 1 and can be obtained using standard techniques; 
we include the analytic expression in appendix B. 

vp -»- •-► 

I 


-<- •- tr 

Figure 1. Leading DM-SM interactions for the case where the effective vertices (represented by 
black circles) are generated by neutral fermionic mediators. The t-channel diagram (<& exchange) 
generates the leading contribution oc Iciiil"*; the s-channel diagram {Z exchange) generates a resonant 
contribution oc that is significant only when ~ mzl‘2.', see (5.5). 



The cross section for Tff' annihilation into heavier fermions has two resonant contributions 
from the Higgs and the Z boson generated by the diagrams in figure 2; their expressions are 
also given in appendix B. From these results, and using the approximations described in [28], 
we readily obtain. 
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where in the first line we ignored 0(m|,/A^) corrections and 
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(5.2) 

(5.3) 

(5.4) 

(5.5) 
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(5.6) 


The expressions (5. 3-5. 5) correspond to the s-wave annihilation processes for the correspond¬ 
ing channels. 

Using standard results [28] we use these expressions to derive the relic abundance: 


where Mpi denotes the Planck mass, 5 * 5 , g'* denote, respectively, the relativistic degrees of 
freedom associated with the entropy and energy density, and 


Xf = = In (0.152m^^) — - In [In (0.152m,i„f)] , 

Tf 2 


and 77 is the freeze-out temperature. 


(5.8) 



Figure 2. Resonant contributions for the —> // annihilation cross section. 

The expression for fl can now be compared to the result inferred from CMB data obtained 
by the Planck experiment [10]: 

f^Pianck/i^ = 0.1198 ± 0.0026 (3cj). (5.9) 

Outside the resonance region Q is determined by the T'k —)• pp cross section (5.5) and so will 
be a function of Aeg and mq,; accordingly (5.9) selects a narrow region in the (m>i., Agg) plane 
(see figure 3), which is well approximated by the relation 

Aes — \ —^TeV; ~ 74GeV (non-resonant region). (5.10) 

V mq, 

In addition to the above analytic calculation, we also derived numerically the constraints 
on the model parameters. This calculation was done by selecting 2 x 10^ points in the 7- 
dimensional parameter space {cn, Cm, A, niq,, m$} within the ranges 

1 GeV < mq, < 199 GeV, 1 TeV < A < 5 TeV , 11 GeV < < 836 GeV, 

0 < c„i < 4, [cii] < 10 , -8 < < 0 , -10 < < 0 . (5.11) 
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Figure 3. Constraints on the fermionic DM model derived from the relic abundance inferred from the 
Plank data (at Scr): the light blue region delimited by the yellow line is the allowed region calculated 
using the MicrOmegas code; the thick blue is obtained using (5.7) with (for illustration purposes) 
Cn = 0.08, Cm = 1.5, = —5, = —5 (the effective scale Aes is dehned in (5.6)). The graph 

clearly exhibits the two resonant peaks at = 771^/2 and ruip = TOij/2. 


while keeping miji < m$. For each point fi was obtained using the public codes MicrOmegas 
[29] and CALCHEP [30] (model implementation for CALCHEP was done using the FeynRules 
package [31]). The results are also presented in figure 3 together with the comparison to the 
analytic expressions. 

5.2 Direct detection 

At present the most stringent limit on spin-independent scattering cross sections of DM- 
nucleon particles comes from the LUX experiment [8]. In order to derive the implications for 
the effective theory under study we obtained the DM-nucleon TAA —)■ TAA cross sections in 
the limit where the relative velocity u —)• 0 (we again use MicrOmegas). In this non-relativistic 
limit the elastic amplitudes are divided into spin-independent interactions (generated by scalar 
and vector couplings) or spin-dependent interactions (generated by axial-vector couplings) . 

The terms in (4.1) that are relevant for direct detection are all contained in 
specifically the terms proportional to Cn and to 

(c'lf'P.+ cK!«>Ph) 5.+... (5.12) 

where we used (5.1); the first and second terms generate DM-SM interactions via Higgs and 
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Figure 4. Shaded region: DM-nucleon cross sections for the parameter ranges (5.11) of the neutrino- 
portal DM scenario, all points satisfy the relic abundance constraint (the discontinuous boundary is 
due to limitations in parameter sampling, except for the Higgs resonance peak at ~ 60GeV). 
Black line: limit form the LUX experiment at 90% C.L. [8]; red line: ATLAS limit [32] derived from 
the invisible decay of the Higgs at 90% C.L. 


Z exchanges, respectively. From this expression we extract 

,,3 ,,2 


en = 


IGvr^Aml^ 


Cii 


ez = - 


jm 

CviI 


MR) 

Cvil 


IGvr^A^ 


(5.13) 


that provide estimates of the strength of the H and Z exchanges to the total direct-detection 
cross section. The final result takes the form 


— /^red I (■^nucl/-^nucl)*^p T (1 ■^nucl/-^nucl)*^ri 


TT 


(5.14) 


.p^n denote, respectively the amplitudes for proton and neutron scattering (in units 


where Al- 

of 1/mass^), Ajiucb ^nud denote the atomic number and nuclear charge respectively, and fi^ed 
the A/" — 'll reduced mass 

The resulting spin-independent cross sections are plotted in Fig. 4 for the parameters 
that satisfy the relic abundance constraint (5.9) and lie within the ranges (5.11). We see 
that there are significant regions in parameter space allowed by the LUX constraint [ 8 ] (the 
ATLAS constraint in the figure is discussed below). In figure 5 we plot the effective couplings 
en, ez defined in (5.13) that satisfy both LUX and PLANCK results. 

For the neutrino portal the spin-dependent cross sections are generated by the axial Z 
couplings in (5.12) and are of the same order as the spin-independent ones. Super-K [ 6 ] 
and ICECUBE [7] have published the strongest limits on this cross section, but these are 
still several orders of magnitude above the ones generated by (5.12) and put no significant 
constraints on the model parameters. 
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Figure 5. Allowed values of Ch and ey by the LUX experiment (at 90% C.L.) and the PLANCK 
constraints from the CMB. 


5.3 Indirect detection 


Dark matter particles in the galactic halo have a finite probability to be elastically scattered 
by a nucleus and become subsequently trapped in the gravitational well of an astronomical 
object such as the Sun or Earth. These DM particles will undergo subsequent scatterings, 
until they thermalize and accumulate at the core of the object [1]. The accumulated DM 
particles in the inner core of the Sun or Earth can annihilate into SM particles that can be 
detected, among others, in astrophysical high energy neutrino experiments. 

In the neutrino-portal scenario neutrinos are the most abundant products of DM self- 
annihilation; the reaction occurs through the exchange of dark scalars into neutrinos (cf. 
fig 1). Given the small DM velocities the neutrino spectrum is essentially a delta function 
centered around ~ the corresponding spectrum detected on ground experiments is 
given by [33-35]: 


dN^ 

dEiy 


dvrii^ 


6{E, 


m>Ir) 


(5.15) 


where is the DM-DM annihilation rate, and R is the distance from the neutrino 

source to the detector (Sun-Earth distance or Earth radius for neutrino annihilation in the 
Sun or Earth, respectively). The annihilation rate is determined by the capture rate 
when capture and annihilation processes reach equilibrium in a time-scale much smaller than 
the age of the body (e.g. Sun or Earth), then ~ C^l2. The capture rate depends 

on the DM-nucleus cross sections, the DM velocity dispersion and the DM local density [35], 
roughly C^, oc where the first factor is generated by (5.12). 

Once produced the neutrinos will sometimes convert into muons as they move through 
the Earth, generating a upward-moving muon flux. We use the MicrOmegas package to esti- 


- 13 - 








Figure 6. Blue shaded region shows the dark matter induced upward muon flux from the inner Solar 
core DM annihilation into neutrinos in the neutrino portal effective model. The red line is the Super-K 
limit over WIMP-induced (neutralino-like) upward muon flux from the Sun. 


mate both the neutrino and the upward muon fluxes. The calculation takes into account not 
only the dominating process (5.15), but all TT —)• SM SM channels, and uses the tabulated 
neutrino spectra functions [36] , taking into account effects induced by oscillation and attenua¬ 
tion processes. Figure 6 show the resulting upward muon flux for the parameter ranges (5.11) 
within the effective model, and, for comparison, the SuperKamiokande [37] WIMP-induced 
(neutralino-like) upward muon flux from inner Solar core dark matter annihilation into SM 
particles. The graph shows that within the neutrino portal scenario this neutrino flux is, in 
general, much smaller than what is expected in generic WIMP scenarios. 

The galactic halo is another source of DM annihilation products that may be accessible 
to indirect detection experiments. While the neutrino flux from DM annihilation captured 
in the Sun or Earth depends on the DM-nucleon cross section, the galactic neutrino flux 
depends on the self-annihilation cross section [38, 39]. ICECUBE experiment measures the 
characteristic anisotropic flux of highly energetic neutrinos for different annihilation channels, 
including direct annihilation into neutrinos [40].; however, as shown in figure 7, the experi¬ 
mental sensitivity is still several orders of magnitude below the predictions of the model. 

5.4 Collider onstraints 

In this section we briefly cover the main existing collider experimental constraints on the 
effective-theory realization of the neutrino-portal scenario. These constraints are derived 
from the invisible decays of the Z and Higgs generated by the reactions, Z —)■ 'kT and 
H —)• fl/T (assuming they are kinematically allowed). 
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Figure 7. Blue shaded region shows the DM annihilation cross section 'P'l' —vv within the neutrino- 
portal scenario. The red line is the ICECUBE limit (published results do not extend to masses below 
lOOGeV) 


Z invisible decay. For mz > 2m^ 
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where AxeV = A/(lTeV). Using (5.11) we find that BTZ{Z —> 'h'h) is at most of 
order 5 x 10“^, which is significantly below the error in the invisible branching ratio 
BTZ{Z —)• inv) = (20 ± 0.06)%. This process does not impose a significant restriction 
over the model parameters. 


H invisible decay. There are two processes that contribute: H —)• 'h'h, (the latter 

followed by the prompt decay —)• 'hi/). 

For the hrst case, after a straightforward calculation we hnd, for mn > 2m^, 
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(5.17) 


where we neglected a small contribution from the class VII operators. This then implies 


BTZ{H 'h'h) ~ 3 X 10” 
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where AxeV = A/(lTeV). The direct search for Higgs invisible decays in the ATLAS 
experiment [32] sets an upper limit of 63% at 90% C.L. This limit on BTZ{H —)■ inv) 
translates into an upper limit on the DM-nucleon spin-independent scattering and is 
shown in figure 4. There are no exclusion regions over the natural ranges of the effective 
couplings. 


For the second case we obtain 


where 
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(5.20) 


The corresponding branching ratio is below 0.4% for m-q, > 5GeV, and it again imposes no 
significant constraint on the model. 


6 Conclusions 

We discussed a general effective-theory description of a multi-component dark sector and its 
interactions with the SM generated by the exchange of heavy mediators. We obtained all 
effective interactions generated by operators of dimension < 6 for the case where all dark 
fields transform non-trivially under an unspecified symmetry group ^dm and are all singlets 
under the SM gauge group ^sMj while all SM fields are assumed to be ^dm singlets. We then 
specialized to the case where the underlying theory is weakly coupled and renormalizable, 
and when the mediators are neutral under ^dm x ^sm and are either scalars or fermions. In 
this case we argued that only a relatively small number of the effective operators are relevant 
given the sensitivity of current and near-future experiments. 

We then consdiered in some detail the neutrino-portal scenario where only fermion me¬ 
diators are present, and where the dark sector consists of fermions and scalars such that the 
lightest dark particle is a fermion. This scenario is characterized by having naturally sup¬ 
pressed couplings of the DM to all SM particles except the neutrinos, and has received little 
attention in the literature. We find that there is a wide region in parameter space allowed by 
the current experimental constraints (relic abundance, direct and indirect detection limits). 
In particular, the DM mass is unconstrained for a significant ranges of the remaining model 
parameters; still, an improvement in one order of magnitude in the experimental sensitivity 
would exclude DM masses below ~ m/f/2. We also considered the possible collider con¬ 
straints on the model parameters and found them to be very weak. The cleanest signature 
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of this scenario is the presence of a monochromatic neutrino line, from both the Sun and the 
halo, with energy equal to that of the DM mass, experimental sensitivity would have to be 
improved significantly before this can be probed. 

Though in this publication we concentrated on the effective theory approach to the inter¬ 
actions between the dark and standard sectors, the neutrino portal scenario is clearly easily 
studied using specific model realizations. We will present such a discussion in an upcoming 
publication. 

A The operators 

In this appendix we list the operators referred to in table 1. Using gauge invariance we 
obtain 

G {|^|^ nl^p, Bl,, (A.i) 

where il: = i, q] if)' = u, d, e; ip = cj), cj) {q, i denote the left-handed quark and lepton SM 
isodoublets respectively; u, d the right-handed quark isosinglets, and e the corresponding 
right-handed lepton isosinglet; cj) = ia24>*■, where (T 2 is the usual Pauli matrix). The operators 
involving fermions correspond to the Yukawa terms in the SM. Operators such as “ipi and 
|D())p were not included because they generate redundant terms in the effective Lagrangian 
after application of the equivalence theorem [41] . 

Similarly 

e (A.2) 

and 

e ^ PL,R'f} (A.3) 

where some terms may be absent for specific choices of the dark symmetry group ^dm- Some 
possibilities such as — (/x o zz) have been eliminated by applying the equivalence 

theorem. 


B Cross sections 


In this appendix we provide, for completeness, the expressions of the various cross sections 
used int he calculation of the relic abundance. 


Neutrino final states 
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where 
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Quark and charged lepton final states 
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